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Pulsar bow-shock nebulae 

III. Inclusion of a neutral component, and Ha luminosity. 
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Abstract. The interaction of the wind from a pulsar (or more generally from a star) with the ambient medium 
gives rise to the formation of a bow-shock nebula. We present a model of adiabatic bow-shock nebulae, including 
the presence of a neutral component in the interstellar medium. As demonstrated in a previous paper (Bucciantini 
& Bandiera 



2001 



hereafter Paper I), the velocity and the luminosity of the pulsar, as well as the characteristics 
of the interstellar medium (ISM), play an essential role in determining the properties of bow shock nebulae. 
In particular, the Balmer emission is strongly affected by the rates of hydrogen charge-exchange and coUisional 
excitation. So far, only one-component models have been proposed, treating the neutrals as a small perturbation. 
The distribution of neutral hydrogen in the nebula, derived by our model, is then used to determine how the Ha 
luminosity varies along the bow shock. We find that the luminosity trend is different for charge-exchange and 
coUisional excitation, and that processes like diffusion or ionization may reduce the emission in the head of the 
nebula. 
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1. Introduction 

Pulsars are known to be sources o f ultr arelativistic mag- 
netized winds (see e.g. Michel & Li 19991) . If during the su- 
pernova explosion the newly formed neutron star acquires 
a high velocity, it may escape from the supernova remnant 
before its pulsar activity has ceased. In such a case the pul- 
sar wind interacts with the interstellar medium (thereafter 
ISM), possibly forming a steady flow, with a bow shock 
in the direction of the stellar motion. So far only four 
such nebulae are known: PSR 1957+20 (Kulkarni & Hester 
198^ ), PSR 2 224+65 (Cordes et al. [l99^, C hatterjee & 



Cordes |2002| ), PSR J 0437+ 4715 (Bell |1995D , and PSR 
0740-28 (Jones et al. ^002[ ). Probably even the nebula 
associated with the neutron star RX J185635-3754 (van 
Kerkwijk & Kulkarni 2001) may be included in the list, 
even if there is no agreement whether the nebula is a bow 
shock or the result of photoionization. 

All the nebulae discovered have been detected in 
Balmer lines (mostly Ha). A pure Balmer spectrum is 
the signature of a non-radiative shock moving through 
a partially neutral medium (Chevalier & Raymond 198C, 
Chatterjee & Cordes 2002). When recombination times 
are long the relevant processes for optical emission are ex- 



citing charge- exchange and coUisional excitation of neutral 
H atoms. 

As mentioned in a previous paper of this series 
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(Bucciantini 2002, hereafter Paper II), the ISM is seen by 
the pulsar as a plane parallel flow with a constant density 
(at least on the length scale of the bow shock). If the ISM 
is partially neutral, and if the neutral atoms have mean 
free paths comparable with, or longer than the size of the 
bow shock a pure fluid treatment is invalid. The interac- 
tion of the relativistic magnetized pulsar wind with the 
ionized component of the ISM is mediated by the mag- 
netic field advected by the pulsar wind and compressed in 
the head of the nebula. The neutral atoms interact with 
the ions (mainly protons) thermalized in the bow shock. 
The penetration thickness for the H atoms changes from 
pulsar to pulsar depending on their velocity, their lumi- 
nosity as well as the density and the neutral fraction of the 
ISM; and the same holds for the diffusion and ionization 
(Paper I). Thus to have a more realistic picture of bow- 
shock nebulae in this regime, it is necessary to include the 
effect of mass loading from a neutral component as well 
as the diffusion and ionization of dragged atoms. 

As a starting point we have used a two thin layers 
bow shock model, for which an an alytic solution may be 
derived (Comeron & Kaper 199S). In Paper II we have 
verified, using a numerical hydrodynamic code, that this 



2 



N. Bucciantini: Pulsar bow-shock nebulae 



analytic model gives reasonable estimates for the values 
of surface density and tangential velocity in the external 
layer. These estimates describe the fluid structure in the 
external layer and have been used to model the interaction 
and the evolution of the neutral atoms. 

Once the density of the neutral H atoms in the exter- 
nal layer has been derived, it is used to evaluate the Ha 
luminosity of the nebula, and its variation along the bow- 
shock. The comparison of the observed luminosity with 
that obtained from a thin layer model may be used to de- 
termine the properties of the local ISM, as well as those 
of the pulsar wind. In Section 2 we present the equations 
we have used to derive the distribution of neutral atoms, 
and the solutions derived integrating these equations; in 
Section 3 we derive the Ha luminosity of the nebula and 
discuss what informations may be obtained from observa- 
tions. 



2. The model 

2.1. The zeroth order two thin layer solution, and the 
effect of dragging charge-exchange 

As said in Paper I, the parameters defining the scenario of 
a runaway pulsar moving supersonically in the ISM are: 
the pulsar energy losses, L; its velocity in the rest frame 
of the ambient medium, Va] the density respectively of the 
ionized and neutral component of the ambient medium, 
rio, SoTio. A supersonic motion implies that the ram pres- 
sure of the ISM is greater than its thermal pressure, so 
that the latter may be neglected in the model. For typical 
pulsar velocities the related Mach number Mq is <L 4 (for 
a diffuse warm ionized medium). 

In Paper II we stressed the fact that a simple bow- 
shock model (Wilkin 1996 ), may not be adequate for the 
case of a pulsar. The pulsar wind is magnetized; thus the 
basic assumption of the Wilkin solution (efficient mixing 
between stellar wind and ambient medium) is not satisfied. 
For a better description we need to follow separately the 
two layers, one containing the shocked stellar wind and the 
other the shocked external medium, allowing two different 
tangential velocities. In this work we are mainly interested 
in the behaviour of the external layer of thermalized ISM. 



If we retain the thin layer assumption, equations de- 
scribing the shape, tangential velocity and surface density, 
can be derived from the momentum and mass conserva- 
tion. We shall start from the freely crossing case (the neu- 
trals do not interact), then adding on it the effect of mass 
loading from a neutral component. In fact in typical cases 
the density of neutral hydrogen may be of the same order 
as that of protons. 

Let us first summarize the freely crossing case. With 
reference to Fig. [|, let us consider a point (P) on the 
layer with distance r from the star, where the incident 
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Fig. 1. Geometry of the two thin layer model (completely 
passing case) . Jo and no are the momentum flux, and the 
numerical density of the ambient medium; and L are 
the stellar wind momentum flux and the pulsar luminos- 
ity; Vo is the pulsar velocity with respect to the external 
medium. 

fluxes of momentum are Jo and J*. The conservation of 
perpendicular momentum gives: 



Jo sin ai — J, sin a2 = 



R 



(1) 



where Ufi is the surface density of the ionized component 
in the external layer, Vf is the tangential velocity in the 
external layer, R is the curvature radius of the front, and 
is the momentum flux in the internal layer. For the 
tangential component of the fluxes in the external and 
internal layer we have respectively: 



Jo sinai cosai — 



J^ sin a2 cos a2 — 



ds ' 
ds ' 



(2) 



(3) 



where ds is the line element. Finally the mass conservation 
equation for the external layer, where surface density is 
required for the evaluation of penetration thickness, reads: 



af^Vf 



2vo 



y- 



(4) 



These equations may be written in a simpler form using 
dimensionless variables. The lengths y and x scale with 
the distance to the stagnation point (Cordes 1996): 



do = \J J*/ Jo = yjLl^-KCnaVl, 



(5) 



while the velocity and surface density in the external layer 
are respectively scaled with Va and no do , and the momen- 
tum flux in the internal layer with Jo do- Equations to ^ 



N. Bucciantini: Pulsar bow-shock nebulae 



3 



may be written: 

dz 

^ dy 



n{yz 



= 0, 



,,21 



= ^' 

ydy n 



(6) 



djy^r) _ {xz + y){yz - x) 
ydy 



0, 
= 0. 



Where z — dx/dy and n = (1 + z^)^^^^ . This set of equa- 
tions can be integrated numericaUy, starting from the stag- 
nation point, and using a Taylor expansion near the apex. 
It is remarkable to see that these equations do not contain 
any free parameter, so there is only one solution in these 
dimensionless coordinates (solid lines of Figs. 2 and 3). 

Once the freely crossing solution has been obtained 
from the integration of Eqs. H, it is possible to introduce 
the effects of a neutral component of ISM (that we assume 
to be mainly composed by neutral hydrogen). As said in 
Paper I, for typical pulsar velocities the main process re- 
sponsible for the dragging of neutral H in the external 
layer is charge-exchange. Consequences of the presence of 
a neutral component are: 1- a new source of perpendicu- 
lar and tangential momentum, 2- the need for a further 
mass conservation equation for the neutrals. We assume 
that, once a neutral atom has been captured in the exter- 
nal layer, it can neither escape via diffusion nor be ionized 
(an assumption whose implications will be analyzed in sec- 
tion 2.2), and that there is no shear between the neutral 
and ionized material. The penetration thickness of the ex- 
ternal layer is not constant but tends to increase moving 
away from the stagnation point (Paper II). The fraction 
of neutrals stopped in this layer is: 



1 - e^p{-acx{Vrel)afiVrel/Vo sinofi). 



where acx is the cross section of charge- exchange, Vrei is 
the relative velocity between the incoming H atoms and 
the hot shocked plasma and we have assumed a constant 
state (i.e. no variations of density, temperature, velocity) 
across the external layer. We want to stress that the con- 
stant state assumption is valid only near the head of the 
bow shock where the geometrical thickness of the layer is 
small compared to the size of the nebula. In order to eval- 
uate the relative velocity we assume that all the ionized 
material in the external layer moves tangentially with a 
speed Vf: this assumption is good (Paper I) if the penetra- 
tion thickness is not too high, otherwise it would be better 
to use the post shock velocity (according to the Rankine- 
Hougoniot jump conditions). However the use of the post 
shock velocity is justified only in case F 1, anyway giv- 
ing results not too different from the other case. The set 
of Eqs. |6[ becomes: 

// ^ -? ^ ^dz „ n(yz — xY 



dy 

d{y{af^+afn)v)) 
ydy 



-(1 + Soi^) = 0, (8) 



djy^r) _ {xz + y){yz - x) 
ydy 



d{yafiVf) 



= 0, 



dy 

diyafnVf) 
dy 



y^o, 



y^oF = 0, 



where a/„ is the surface density of the neutral component 
in the external layer. It is convenient to parametrize F 
using a power law for the cross section in the typical range 
of pulsar velocities (Newman & al. 1982, McClure 1966, 
Peek|l966D: 



Vrel 
Vrel.O 



-0.37 



(where we have taken Vrei,o — 1) so that: 



F=l-F^ 



Ke,(j,)«-63<T^,n/1.25) 



(9) 



(10) 



(7) S, 



Therefor 1 — Fq corresponds to the fraction of neutrals 
captured at the stagnation point, in the fully ionized case, 
which defines the properties of the true physical nebula 
(the factor 1.25 corresponds to the surface density of the 
external layer on the axis, in the freely crossing case). 

Contrary to the freely crossing case, Eqs. ^ cannot be 
integrated from the head, because they are not closed (i.e. 
in the Taylor expansion of the equations the terms of order 
0(2/") are function of those of order 0(1/"+^)). Instead we 
have used a perturbative method: from the freely cross- 
ing solution, we evaluate a first approximation of F; this 
approximation is then used in Eqs. ^ instead of the self- 
consistent expression, both to obtain the Taylor expansion 
and to integrate numerically the equations; we repeat it- 
eratively till F converges. 

In Fig. 2 we present the solution of Eqs. || for the case 
= 1 (an intermediate case between the highly ionized 
and highly neutral ISM) for different penetration thick- 
nesses in the head. The solid line represents the freely 
crossing case (it is a good approximation if the thickness 
remains ^ 1 all along the nebula), while the dotted line 
refers to the completely dragging case (which is simply 
a bow shock scaled for a greater external density). The 
stagnation point distance from the star decreases with in- 
creasing values of Fo, and its value is: 



d ~ do/V(l + 2o(l--Fo)). 



(11) 



The self-consistency level of our solution decreases moving 
to the tail of the nebula, away from the star, because the 
thin layer approximation is no longer acceptable, an so 
the constant (perpendicularly to the bow-shock) internal 
state assumption lose its validity. 

Fig. 3 shows the surface density of the ionized compo- 
nent in theSo = 1 case: also here the solid line represents 
the freely crossing case while the dotted line is the com- 
pletely dragging solution. It is evident that the solution 
with highest F tends to the dotted line in the tail of the 
nebula. Fig. 4 gives the surface density of the neutral com- 
ponent in the various cases obtained from Eqs. 0. Once the 
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Fig. 2. Shape of the bow shock for various penetration 
thickness in the head. The sohd Hne is the freely crossing 
case; the dotted one is the completely dragging case. The 
short-dashed, long-dashed and dot-dashed lines refer re- 
spectively to three values (0.1, 0.3, 1.) of the parameter 
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Fig. 3. Surface density of the ionized component in the 
external layer the bow shock for various penetration thick- 
ness in the head (So = 1). The various lines refer to the 
cases of Fig. 2. 



neutral and the ionized component density has been ob- 
tained it is quite easy to derive the Ha luminosity along 
the nebula. However, before doing this it is necessary to 
consider the effect of other processes like diffusion and 
ionization (see Paper I) which modify the neutral distri- 
bution. 
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Fig. 4. Surface density of the neutral component in the 
external layer the the bow shock for various values of the 
penetration thickness in the head (So = 1). The various 
lines refer to the cases of Fig. 2. 

2.2. Importance of diffusion and ionization and their 
inclusion in the model 

As previously mentioned the neutral H atoms, once 
stopped by charge- exchange with thermalized protons of 
the external layer, may diffuse or be ionized as a conse- 
quence of processes of charge-exchange and ionization in 
the internal hot plasma. Before introducing these effects 
in the equations it is necessary to discuss the physical 
regimes of the pulsar bow shock nebula, and their varia- 
tions along the nebula. Comparing the time scale of the 
flow Tfiux ~ do/vo with, respectively, the time scale of 
diffusion and ionization, tells what are the importance of 
these two processes in the various pulsar cases. 
The time scale for diffusion is: 

A cTfiacxivth) 



Vth 



(12) 



where A is the geometrical thickness of the external 
layer, the cross section of diffusion acx is that of charge- 
exchange, and the thermal velocity vth is that obtained 
from the Rankine-Hugoniot condition on the bow shock 
^ 0.8 Wo sinai. Moving away from the head this time scale 
tends to increase, so that if it is greater than the flow time 
on the head, diffusion may not be taken into account. The 
time scale for ionization is: 

Tion ^ '/ ^ , (13) 

[VthjVth 

where (T,;o„ is the cross section of ionization, that, in the 
range of pulsar velocities, can be approximated by the 



relation (Ptak & Stoner |T97|) : 

/ Vth 



i{vth) 



i{vth,o) 



\Vth,0 



2.9 



(14) 
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Also this time scale increases in the tail of the nebula, 
mainly as a consequence of the decreasing of the thermal 
velocity. 

In the case of very luminous pulsars it is also neces- 
sary to determine the relative importance of diffusion and 
ionization, because atoms are ionized before the may dif- 
fuse (Paper I) . The ratio r^i / / /rion tends to decrease away 
from the head because the cross section of ionization has a 
steeper velocity dependence than that of charge-exchange. 

The effect of diffusion and ionization may be included 
in the equations to give a more realistic model, but some 
simplifications are necessary. We have assumed that atoms 
in the external layer have the same probability to diffuse 
in the internal stellar wind region as to diffuse in the exter- 
nal region of unperturbed ISM. In reality there may be an 
asymmetry in the diffusion not taken into account in our 
model, however the results of its inclusion are not signi- 
ficatively different, both in terms of surface densities and 
shape of the nebula. Diffusion acts in two ways: 1- reduc- 
ing the surface density of neutral atoms, 2- bringing away 
from the layer an amount of tangential momentum, pro- 
portional to the number of diffused atoms. To treat prop- 
erly diffusion we need to know the distribution of neutrals 
in the layer, but thin-layer models lack this piece of infor- 
mation. Contrary to diffusion, ionization does not change 
the total surface density in the external layer, but simply 
changes neutrals into ions. This may modify the geometry 
of the nebula because the increased ionized surface den- 
sity makes the layer thicker to the penetration of neutrals. 
Using Eqs. ^ and |lj for the parametrization of the cross 
sections we find that the rate of diffusion C and ionization 
B may be written: 



2 r 



C = Co- 



B = So ^ 



(15) 
(16) 



where all the unknown parameters are enclosed in the 
quantities Co and Bo, and Ao is the thickness of the 
layer along the axis in the freely crossing case. From a 
numerical simulation (Paper II) we find the parametriza- 
tion A = Aod(l + y^/5)/do. Eqs. § then become: 

{{aj, + aj^)v} + $,)- - n(l + E^F) + ' = 0, 

d(,(a,.+a,>P _. 
ydy n 



djy^r) _ jxz + y){yz ~ x) 
y dy nr^ 



= 0, 



dijjUfiVf) 



y{l + B)^Q, 



djyCTfnVf) 

dy 



dy 

- y(So F-C-B) 



0. 



They can be integrated in the same way as discussed in 
the previous subsection, in fact the quantities C and B 
may be evaluated self-consistently. 
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Fig. 5. The effect of diffusion on the surface density for the 
neutral component {Xio = 1). Panels refer to: A to com- 
pletely dragging case, B C and D respectively to the values 
1, 0.3, 0.1 of the parameter acx{vrei.o)nodo- Continuous 
lines are the solutions without diffusion, dotted lines have 
Co = 0.25, dashed lines have Co = 0.8, dash-dotted lines 
have Co = 4. 

Figs. 5 and 6 show respectively the effect of diffusion 
and of ionization on the surface density of the neutral 
component in the external layer, for various values of the 
coefficient Co and B^. Both diffusion and ionization tend 
to reduce the density of H atoms. However this effect is re- 
duced in the tail: diffusion decreases both as a consequence 
of the greater thickness of the layer and as a conseguence 
of the reduction of the mean free path of neutral H atoms; 
ionization simply vanishes because the cross section de- 
creases rapidly with the thermal velocity. Fig. 7 shows the 
surface density of the ionized component: in case of hight 
ionization rates it tend to decrease in the head with re- 
spect to the values on the axis. This influences the Ha 
luminosity, making it less that what may be derived from 
a simple case with only dragging charge-exchange. 

3. Ha luminosity 

Once a neutral component has been introduced in the 
model, it may be used to evaluate the Ha luminosity of the 
bow shock nebula and its variation moving away from the 
head. In Paper I, we have given an approximate estimate 
of such a luminosity, using some simplifications whose va- 
lidity were limited to the stagnation point. The knowledge 
of the global distribution of neutrals and ions in the nebula 
makes it possible to evaluate the variation of luminosity 
along the bow shock. 
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Fig. 6. The effect of ionization on the surface density for 
the neutral component (Sq = 1). Panels refer to: A to 
completely dragging case, B C and D respectively to the 
values 1, 0.3, 0.1 of the parameter acx{vrei,o)'nodo- Solid 
lines are the solutions without ionization, dotted lines have 
-Bo = 4, dashed lines have Bo = 1, dash-dotted lines have 
Bo = 0.25. 




Fig. 7. The effect of ionization on the surface density for 
the ionized component (So = 1). Panels refer to: A to 
completely dragging case, B C and D respectively to the 
values 1, 0.3, 0.1 of the parameter crcx{vreifl)n'odo- Solid 
lines are the solutions without ionization, dotted lines have 
Bo = 4:, dashed lines have Bo = 1, dash- dotted lines have 
Bo = 0.25. 



Before investigating the effect of diffusion and charge- 
exchange we need to discuss the distribution of neutrals 
in the external layer. In those cases in which the pene- 
tration thickness for H atoms if not too high, we assume 
that the neutrals are homogeneously distributed in the 
external layer so that the luminosity per unit volume is 
constant in the direction normal to the layer. But in the 
highly dragging case it is possible that neutrals cannot 
reach the contact discontinuity, and remain confined in a 
thinner shell near the bow shock (al least in the head of 
the nebula) . Anyway only a very high resolution mapping 
may discriminate among the above cases. In fact the in- 
strumental resolution may spread the pattern. This is a 
problem if the signal to noise ratio is low so that only the 
more luminous regions of the nebula are detected. 

To investigate properly the emission from such neb- 
ulae we have to consider the various processes responsi- 
ble of the excitation of H atoms. They may be divided in 
two types: primary processes, which involve neutral atoms 
coming directly from the ambient medium; secondary pro- 
cesses, related to atoms dragged and thermalized in the 
external layer. In a non-radiative shock moving through a 
partially neutral medium the dominant processes for op- 
tical emission are exciting charge-exchange and collisional 
excitation of neutral H atoms, which give the observed 
Balmer spectrum. However primary processes usually have 
a very low efficiency for incoming neutrals (Paper I) so. 



the nebula can easily be detected only if the emission is 
dominated by the secondaries, otherwise a large amount 
of neutrals in the ISM is required to give enough flux. In 
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Fig. 8. Peak luminosity along the nebula from secondary 

charge-exchange, for the various solution shown in panel 
A Fig. 5. The flux is in unit of geff noO'cx{vrei,o)vodo The 
solid line refers to the completely dragging case without 
diffusion or ionization; dotted dashed and dot-dashed lines 
refer respectively to the case Co = 0.25, 0.8, 4 
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creases in the tail because the rate of interaction is a steep 
function of the temperature, and vanishes in the tail. So 
it is possible to discriminate between a flux dominated 
by charge-exchange or by coUisional processes, analyzing 
the Ha luminosity along the nebula. However even if the 
efficiency of coUisional excitation is one-two orders of mag- 
nitude higher than that of charge- exchange, its cross sec- 
tion is much less at typical pulsar velocities, so that only 
in the faster cases ionization may give a non neglectable 
contribution, at least in the head of the nebula. 
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Fig. 9. Peak luminosity along the nebula from secondary 
coUisional excitation, for the various solution shown in 
panel A Fig.s 6 7. The flux is in unit of euoBoVo- The 
solid line refers to the emission in the completely drag- 
ging case with an ionization coefhcient Bo — A the dashed 
line refers to the case Bq — \, and dash-dotted line to the 
case Bo = 0.25 . 



what follows only secondary charge- exchange and ioniza- 
tion are considered. 

Fig. 1^ shows the behaviour of peak luminosity along 
the nebula for secondary charge- exchange. The various 
curves refer to the same cases of panel A in Fig. 5. There 
are various effects that contribute to the shape of the 
Ha luminosity. The solid line represents an upper limit 
while the others show the effect of diffusion. The luminos- 
ity tends to increase as a consequence of the increasing 
column length of integration, but this effect is reduced 
by the diminishing of the ionized and neutral density, and 
this gives the behaviour of the various lines. While the case 
without diffusion gives a quite constant emission along the 
nebula, the presence of diffusion reduces the Ha luminos- 
ity in the head (where the efficiency of diffusion is grater) , 
and this effect could be used as a tool to determine the 
physical conditions of such nebulae. 

We have not taken into account the presence of a 
threshold for the excitation of H atoms: we considered 
a pulsar so fast with respect to the ambient medium that 
even in the tail the thermal velocity is above the excita- 
tion energy (~ 10 eV). In the case of a slower pulsar, the 
temperature in the tail may not be enough to allow excit- 
ing charge-exchange, so the Ha flux may be lower. This 
is the case of PSR J0437-F4715 (Bell & al. |1995D : only 
the head of the bow shock nebula have been detected so 
far despite the fact that this is the nearest object. In fact 
its velocity is estimated to be ~ 80 km/s, very near the 
threshold, while the other three nebulae, associated with 
faster pulsars, are detected even in the tail. 

Fig. 1^ shows the behaviour of peak luminosity along 
the nebula for secondary coUisional excitation. The various 
lines refer to the same cases of panel A Fig. 5. Only models 
with ionization are shown. The Ha emission rapidly de- 



4. Conclusion 

We have modified the two thin layer solution (Comeron 
& Kaper 1998) for bow shock produced by runaway stars 
(with the aim of deriving a model in the case of pulsar 
wind bow shock nebulae) with the inclusion of a neutral 
component. We find that the ISM neutrals lead only to 
minor effects to the shape of the nebulae (see Fig. 2), 
but plays an essential role in their luminosity. We have 
also introduced in our model the possibility that the H 
atoms, once stopped by the thermalized protons, may dif- 
fuse out of the layer or be ionized, as a consequence of 
charge-exchange and coUisional ionization. Both processes 
act depleting the neutral surface density in the external 
layer, but their efficiency is higher in the head of the neb- 
ula both as a consequence of the decrease of the rate of 
charge-exchange and ionization in the tail, and the in- 
crease of the geometrical thickness of the layer. 

The distributions obtained have been used to evaluate 
the Balmer emission along the nebula. The study of such 
emission can be a good tool to understand the proper- 
ties of local ISM as well as those of the pulsar. Even if 
this is not a simple task, and requires very good measure- 
ments, which only for two objects are actually available, 
an accurate study may show which hypothesis can be ap- 
plied (like, for example, the absence of preionization both 
from the pulsar surface radiation, or from the emission of 
shocked ISM). In Paper I, we stressed the point that the 
total Ha flux depends on the velocity of the pulsar. 

Here we evaluate its variations along the nebula. The 
behaviour is different if the dominant process is exciting 
charge-exchange or coUisional excitation, so not only the 
total flux but also its variation moving away from the stag- 
nation point can be used. Concerning the emission we flnd 
that: 1- charge-excange and coUisional excitation have dif- 
ferent trends; 2- processes like diffusion or ionization may 
play an essential role changing the emission in the head 
of the nebula. This not only allows one to put stronger 
constraints on the regime of such nebulae, but it may be 
used to verify the validity of the model, by comparison 
with the observed cases. 
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